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ABSTRACT 

We present results from a fifteen-month campaign of high-cadence (^3 days) mid-infrared Spitzer 
and optical {B and V) monitoring of the Seyfert 1 galaxy NGG 6418, with the objective of determining 
the characteristic size of the dusty torus in this active galactic nucleus (AGN). We find that the 3.6 /tm 
and 4.5 fdia flux variations lag behind those of the optical continuum by 37.212 2 days and dT.llg/ 
days, respectively. We report a cross-correlation time lag between the 4.5 /rm and 3.6 /im flux of 
13.91 q 1 days. The lags indicate that the dust emitting at 3.6 /rm and 4.5 fim is located at a distance 
1 light-month (« 0.03 pc) from the source of the AGN UV-optical continuum. The reverberation 
radii are consistent with the inferred lower limit to the sublimation radius for pure graphite grains at 
1800 K, but smaller by a factor of ^ 2 than the corresponding lower limit for silicate grains; this is 
similar to what has been found for near-infrared (K-band) lags in other AGN. The 3.6 and 4.5 /im 
reverberation radii fall above the K-band r oc size-luminosity relationship by factors < 2.7 and 
< 3.4, respectively, while the 4.5 /im reverberation radius is only 27% larger than the 3.6 /im radius. 

This is broadly consistent with clumpy torus models, in which individual optically thick clouds emit 
strongly over a broad wavelength range. 

Subject headings: galaxies: active — galaxies: individual (NGG 6418) — galaxies: nuclei — galaxies: 

Seyfert 


1. INTRODUCTION 

In the AGN unification paradigm, direct observation 
of the nucleus is blocked by a toroidal structure of 
dusty molecular gas for a range of viewing angles (e.g., 
lAntonuc^ I1993D . As this dust absorbs UV-optical ra¬ 
diation from the accretion disk and re-emits in the in¬ 
frared (IR), this structure is also thought to be the dom¬ 
inant source of IR radiation in most AGN. Understanding 
this obscuration of the central engine is therefore impor¬ 
tant to understanding the physical processes operating in 
AGN and more generally, their role i n galaxy ev ol ution. 

The ob s ervati on al evidence (lAntonuccil 119931 
iJaffe et ^ 12004 iTristram et al.l I2007D . indicates 
that the obscuring structure is geometrically and 
optically thick, although a warped thin disk that 


extends throughout the ho st galaxy has also been 
proDOsed ljSanders et ahlflo^ . The conventional picture 
is that of a compact, but geometrically thick, torus of 
optically thick molecular clouds with a size of a few 
parsecs (lAntonucci fc Mill^ 119851 iKrolik fc Begelm^ 
119881 iPier fc Kroliklll992D . Models in which the vertical 
thickness is supported by large ra ndom velocities due 
to el astic collisions between clouds (|Krolik fc Begelmaii 
1988), or by IR radiation pressure (jPier fc Krolik 
1992 : iKrolikI 120071). or bv turbule n ce induced by su- 
pemovae (jWada fc Normanl 120021 iSchartmann et all 
1200^ have been explored. In an alternative class of 
models, the dusty material is not part of an essentially 
static torus, but is rather embedded in an outflow¬ 
ing hydromagnetic wind launched from the accretion 
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Table 1 

Observations 


Telescope 

start date 

end date 

# obs 

instrument 

filter 

aperture 

Liverpool Telescope 

08-06-2011 

10-21-2012 

64 

RATCam 

Bessel B 

1.5” 

Faulkes Telescope North 

08-10-2011 

09-30-2012 

60 

fs02 

Bessel B 

1.2” 

SU Fountainwood 0.4-m 

05-19-2012 

12-18-2012 

48 

SBIG ST-8300 

Johnson-Cousins B/V 

3.5” 

Spitzer 

08-01-2011 

01-04-2013 

170 

IRAC 

chi/ch2 

1.8” 


disk i e.g.. iBlandford fc Pavnd Il982t lEmmering et alJ 
— /l2 ■ 


19921: iBo ttorff et alJ 

Dorodnitsvn et alJl2012n . 


19971: lElitzur &: Shlosmani 120061: 


for the torus 
energy distribu- 


Dust radiative transfer models 
broadly reproduce the IR spectral 
tion (SED) of AGN. Of necessity, early radiative 
transfer mo dels assumed smoot h density distribu¬ 
tions ( e.g. [Pier fc KrolikI Il993l: iGraiiato fc Danesd 
11994 lEfstathioii fc Rowan-RobinsonI 1199511 . but 
more recently, models for clu mpy dust dis tribu¬ 
tions have been developed (e.g., iNenkova et al.l 
Dullemond &: van Bem niel 120051: iHonig e t al.l 


2002 


2006 


Schartmann et al.n2008r^enkova et al.ll2008allbl i. These 
“clumpy torus” models are more successful in reproduc¬ 
ing certain details of the SED such as, fo r example, the 
strength of the 10 nm sil icate feature (jNikutta et al.l 
l2009l : 1^nkova et al.l(2008fl l. 

The torus is too small to be directly imaged by 
any existing single telescope. Some constraints on its 
size and structure can be inferred fro m SED-fitting 
using r adiative tr a nsfer m odels (e.g.. INenkova et"^ 
2008bl: i Mor et aP 120091: iHonig fc Kishimotol 2M 


Ramos Almeida et al.l 20111 : lAlonso-Herrero et al.ll20lTl i 


but there are many theoretical and observational un¬ 
certainties which obfuscate the results. Other methods 
are therefore required, the two most important being 
reverberation mapping and, for relatively close objects, 
IR interferometry. 


Fol lowing the seminal work of (|Blandford fc McKeel 
1198211 . the reverberation mapping technique has been 
well developed and extensively applied to studies of the 
broad emission line region (BLR). Time series analysis 
of the response of the broad emission lines to variations 
in the UV or optical continuum (as proxies for the AGN 
ionizing continuum) has revealed the characteristic size 
of the BLR in about 50 AGN, enabl ing estimates o f black 


hole masses and Eddington ratios (lPetersonlll993 


2006 


iGaskelll 120091 : iGalianni fc Horn^ 1201,4 iDu et al.l 12014 . 
and references therein). It has also been determined 
that the BLR follows a size-luminosity relationship of 
the form R oc (iPeterson et ahl 12004 iGreene et al.l 
[Mol: iBentz eT^I^Om ^ 

Near-IR (K-band) versus optical (V-band) reverber¬ 
ation lags have been measured for a ro und 20 Seyfert 
galax ie s (lOknva nskii fc Hornd 120011: iMinezaki et al.1 
2004; iSuganuma et al.l 120061 : iKoshida et al.l 120091 
2014f l. As dust grains emitting in the K-band have 
temperatures close to the sublimation temperature 
1200 — 1800K, depending on grain composition), 
these lags are thought to represent the inner radius of 
the torus. The K-band reverberation lags are found 
to be larger than those of the BLR, while following a 
similar R oc size-luminosity relation, implying that 


the BLR is bounded by the dust distribution, consis¬ 
tent with the central idea of the AGN unification scheme. 

The inner regions of several bright, nearby 
Seyfert galaxies have been direc tly studied using 

near-lR (K-ban d) _ in terferometry (ISwain et al.l 120031; 

Kishimoto et al.ll2009l: IPqtt et al.ll2010HKishimoto et al l 


2014 iWeigelt et al.ll2012h .~ The effective ring radii de¬ 

rived from the observed visibilities scale approximately as 
and are comparable with or slightl y larger than the 
radii derived from reverb e ration lags ( Kishimoto et all 
1201 If) . Since iJaffe et all (120041 ’s pioneering study of 
the archetypal Seyfert 2 galaxy, NGC1068, mid-lR 
(8 — 12/im) interferometric observati ons have also 
been obtained for ft! 20 A G N (e.g.. iTristram et al.1 
2007t iB urtscher et al.l 120091: iKishimoto et al.l 120091 : 


Tristram et al.ll2009l: iHbnig etjlj 2Q131) In a r ecent anal¬ 
ysis of the available data. IBurtscher et al.l (j201,4l find 
that while the mid-IR source size scales with luminosity 
in a manner similar to that seen in the near-IR, the 
inferred size is more than an order of magnitude larger 
than the measured K-band size and the scatter is quite 
large. 

Here we report initial results from a mid-IR (3.6 /rm 
and 4.5 fim) reverberation-mapping campaign using the 
Spitzer Space Telescope in its “warm mission”. Our mo¬ 
tivation is to probe the dust distribution at spatial scales 
intermediate between the innermost regions probed by 
the K-band observations and the outer, cooler regions 
probed by mid-IR interferometry. Furthermore, vari¬ 
ability at 3.6 fj,m and 4.5 fj,m should be less susceptible 
than the 2.2 /rm K-band to complicating effects such as 
dust sublimation (IMinezaki et al.ir2004 IKishimoto et al.l 
[2?fn or contamination by variable accretion d isk emis¬ 
sion (jTomita et al.ll2n0^ IKishimoto et al.ll2007f) . During 
a 2-year campaign, we monitored a sample of 12 Seyfert 
1 AGNs at cadences of 3 and 30 days during the first and 
second year, respectively. We selected our targets based 
on their proximity (z < 0.4) and their location near one 
of Spitzer’s continuous viewing zones. We obtained B 
and V images of the targets over the same period us¬ 
ing the Liverpool Telescope, the Faulkes Telescope North 
and the Southwestern University 0.4-m telescope. 

In this work we describe our analysis of the first 17 
months of measurements o f the Seyfert 1 NGC 6418 
( Veron-Cettv fc yeronll200^ . a Hubble type Sab galaxy 
( Nair fc AbrahamI 120101 ) with an ap parent magnitud e 
g = 14.87 at a redshift of z = 0.0285 ()Ahn et al.ll2014fl . 
It is classified spectroscopically as a Seyfert 1 on the basis 
of a strong, broad Ha emission line, bu t it is otherwise 
domi nated by the stellar continuum (see lRemillard et ^ 
I1993L who described it as an “embedded” AGN). Never¬ 
theless, it is also an X-ray sou rce with a 0.1-2.4 keV lumi- 
nosity of Lx = 10^^-^® erg/s (jAnderson et al.l[2007t) . We 
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selected NGC 6418 out of our sample due to its larger 
than average variations in the Spitzer channels for the 
first year of data; the result of the analysis of the other 
targets will be presented in a future publication. 

We present our observations and describe our methods 
for measuring the light curves in Section [2j In Section |3] 
we describe the time series analysis technique that was 
used to extract the time lags between the 3.6 /rm, 4.5 p,m 
and optical light curves. We discuss the implications of 
our results in Section lU and present our conclusions in 
Section m Details of our photometric measurements and 
a comparison of two methods for determining time lags 
can be found in the appendices. 

2. OBSERVATIONS 

We will discuss the mid-infrared and optical observa¬ 
tions separately. See appendix|A]for a detailed discussion 
of our photometric analysis. 

2.1. Mid-Infrared 

We monitored 12 AGN using the Infrared Array Cam¬ 
era (IRAC) aboard the Spitzer Space Telescope for a pe¬ 
riod of approximately 2 years during Cycles 8 (program 
80120) and 9 (program 90209) of the “warm” mission. 
All objects were observed in both IRAC Channel 1 (3.6 
/rm) and Channel 2 (4.5 ^m). During Cycle 8, repeated 
observations of each object were obtained at intervals of 
3 days. In Cycle 9, a longer cadence was used, with 30 
day intervals between observations. Here we report re¬ 
sults from the Cycle 8 high-cadence monitoring of NGC 
6418. Images of this object were obtained every 3 days 
from 2011 Aug to 2013 Jan, except for a 30-day gap 
in 2011 Dec. Each image had an exposure of 10 sec¬ 
onds. All the resulting IRAC im ages were mosaiced us¬ 
ing MOPEX (jJacob et al.ll200^ directly from the Basic 
Calibrated Data (BCD) level 1 products. Photometry 
was extracted from the BCD mosaics generated by the 
MOPEX standard pipeline, as described in Section lOl 

2.2. Optical 

Contemporaneous optical monitoring was performed in 
the B and V bands with three ground-based telescopes: 
Bessel B images were obtained with the 2-m Liverpool 
Telescope (LT) on La Palma and the 2-m Paulkes Tele¬ 
scope North (FTN) on Maui; Johnson-Cousins B and 
V images were obtained with the 0.4-m telescope at 
Southwestern University’s (SU) Fountainwood Observa¬ 
tory (see Tabled]). It was not possible to coordinate these 
observations with each other or with the Spitzer obser¬ 
vations, but together they approximately span the time 
period covered by the Spitzer campaign except during 
November 2011, when NGC 6418 was unobservable from 
the ground. The start and end dates of the observations 
with each telescope are given in Table dl 

The exposure times for the optical observations range 
from 60 to 180 sec. Dark/bias subtraction and flat- 
field division of all images from S U were performed us¬ 
ing th e XVISTA software package (jTreffers fc Richmondl 
Il989t) . Images from the RATCam instrument at the LT 
were bias ed subtracted and flat fielded by an automatic 
pipeline (iSteele et al.l 1200411 , as were images taken by 
the FTN. When more than a single exposure per night 
was available from LT and FTN, we stacked and regis¬ 
tered the images using MATCH, an implementation of 


Table 2 

Mean Flux Density 


name 

mean flux density 

3.6 //m 

3.62 mjy 

4.5 /i.m 

3.54 mJy 

SU B band 

0.53 mjy 

SU V band 

1.48 mjy 

LT B band 

0.50 mJy 

FTN B band 

0.54 mJy 


the star matching algorithm of iTaburl (1200711 . and the 
XVISTA package. We then extracted photometry from 
the stacked image. 

Hereafter, we refer to the light curve compiled from 
the LT and FTN observatories as the combined optical 
light curve. The SU observations are used to determine 
the AGN/Host ratio. The mean flux densities measured 
within the aperture in tabled) for all bands are tabulated 
in tabled These flux densities are not host subtracted. 

2.3. Photometry 

The photometric analysis proceeds in two stages for 
the SU dataset : in the first we measure instrumental 
magnitudes for each object (the target plus comparison 
stars) in all exposures; in the second the measurements 
from all exposures in a given passband are combined and 
the measured instrumental magnitude s are subjected t o 
inhomogeneous ensemble photometry (lHonevcuttlll993l. 
For a detailed discussion of these steps see appendix M 
The LT a nd FTN dat asets are reduced using image dif¬ 
ferencing (IAlardll200Clll and references therein. The com¬ 
bined optical and Spitzer light curves are shown in flux 
density, normalized to the mean, in figures [T] and O The 
light curves are also shown after applying a shift equal 
to the time lag computed by the cross-correlation anal¬ 
ysis (Sec. [3]). In figures [T] and [2] the time lag shifts are 
37.2 and 47.1 days for the Spitzer’s 3.6/rm/optical and 
Spitzer’s 4.5/rm/optical, respectively. 

The optical and infrared curves all show clear varia¬ 
tions with similar features on timescales of ~ 100 days, 
but with the variations in the infrared lagging behind 
those in the optical. 

3. TIME SERIES ANALYSIS 

The reverberation lag, t, between the driving optical 
continuum variations and those of the responding IR 
emission gives the characteristic size of the IR emitting 
region. The lag can be determined by cross-correlating 
the two light curves. The application of this technique 
to the broad emission line variability of AGN (“rever¬ 
beration mapping”) is well developed (iGaskell fc Soarkd 
19861: iGaskell fc PetersonlllOSTt lEdelson fc Kroliklll988t 

Maoz fc Netzedll989HKoratkar fc Gaskelllll991h and has 
been widely used to measure the size of the broad line 
region (e.g., IPeterson et al.ll200^ see IPetersonll200l] for 
a tutorial). As already noted, it has also been applied to 
optical and AT-band light curves in order to determine 
the inner radius of the torus (iSugyuma et al.l l2006t 
lOknvanskii et al.ll2006l : iKoshida et al.ll2009ll . 
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Days after MJD 55900 (12-05-2011) 



Figure 1. Spitzer 3.6 p.m and the combined B band optical data 
light curves. The error bars of the 3.6 pm and the combined B band 
optical light curves are the uncertainties reported by MOPEX and 
the image differencing solution, respectively. The bottom panel 
shows the combined optical light curve shifted by +37.2 days. 

We performed cross-correlation analyses for three pairs 
of data sets: 3.6 /im versus combined optical, 4.5 /im 
versus combined optical, and 4.5 /rm versus 3.6 /im. The 
time series analysis was performed between the dates of 
MJD 55900 (12-05-2011) and MJD 56300 (1-08-2013). 
This time span was selected due to the significant optical 
and IR variations of the light curves and because there 
were no large gaps in coverage. For a comprehensive 
and detailed analysis of individual datasets see appendix 
m For each pair, the cross-correlation function (CCF) 
was computed using a lag step size of 1 day. The 
optical observations were not synchronized with the 
Spitzer observations and are typically separated by 
irregular intervals. On the other hand, the Spitzer light 
curves are for the most part more evenly and densely 
sampled than the optical measurements. Therefore, in 
order to compute the IR~optical CCFs, we generate IR 
data points corresponding to the optical observations 




Figure 2. Spitzer 4.5 fim and the combined B band optical data 
light curves. The error bars of the 4.5 fim and the combined B band 
optical light curves are the uncertainties reported by MOPEX and 
the image differencing solution, respectively. The bottom panel 
shows the combined optical light curve shifted by +47.1 days. 


by interpolating within the Spitzer light curves. For 
examples of the CCFs computed for the 3 pairs of light 
curves see appendix [Cl 

The maximum of the CCF yields the lag, t, between 
the two light curves. However, the maximum is not al¬ 
ways well defined, since computed CCFs typically ex¬ 
hibit a broad peak (see appendix [C]) and structure in the 
wings (at large positive or negative lags), which can in¬ 
fluence the calculation of the centroid or mean. A com¬ 
mon approach is to calculate the centroid of the CCF 
using a subset of points whose correlation coefficients 
exceed a certain valu e: for example, 80% of the peak 
value (|Petersonll2001fl . Here, we use a different method 
in which we fit a cubic spline to the CCF and use it to 
set a threshold for the minimum correlation coefficient. 
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Table 3 

Comparison of cross-correlation methods 


name 

3.6 //m-Optical 

4.5 /im-Optical 

3.6 /i.m-4.5 /im 


(lag(day) ± <5) 

(lag(day) ± 5) 

(lag(day) ± <5) 

Peterson et al. 

36.7±3.4 

48.6±3.7 

14.6±6.0 

Zu et al. 

40At°l 

49.5l4,? 

13.211® 

Vazquez et al. 


47.113;] 

13.9l°l 


This minimum correlation coefficient is defined as: 


= max(CCi^(T)) - 2a{CCFf,t{T) - CCF{t)) 

where CCmin is the minimum correlation coeffi¬ 
cient, CCF{t) is the cross-correlation function and 
a{CCFfit{T) — CCF[t)) is the standard deviation of the 
difference between the fitted and actual CCF value. The 
CCF centroid is computed using only values exceeding 
CCmin- For more details see appendix [Cl 

To estimate the uncertainty on the CCF lags, we 
used the cross-correlation centroid distribution (CCCDl 
method (ICaskell fc Peters^ll987HMaoz fc Netzedll989t 

IPeterson et al.1 flOMl . generating 1000 random realiza¬ 

tions of the light curves. The CCCDs for the 3 pairs of 
data sets are shown in Figure [3] and the derived lags are 
listed in Table [3] and in appendix |B] The lag is taken to 
be the median of the distribution and the uncertainty 
is given by the interquartile range. The CCCDs for the 
3.6 p,m versus optical, 4.5 /rm versus optical and 3.6 p,m 
versus 4.5 ^m light curves yield lags of 37.212^2 days 
(31.21^;^ X 10-3 pc), 47.lt|) days (39.5t|^ x lO'^ pc), 
and 13.9t3;8 days (11.71 q; 3 ^ pc), respectively. 

For comparison, we also analyzed our data following 
th e slightly differe n t cros s-corr elation m e thods described 
by IPeterson et al.l (1200411 and iZu et al.l (1201111 . The re¬ 
sults are compared in Table |31 We find that all methods 
yield results which are consistent within the uncertainties 
for all three pairs of light curves. 

4. DISCUSSION 

The dusty torus absorbs UV/optical radiation from 
the accretion disk and re-emits it as infr ared radiation 
(|Telesco et al.l 11984 iSanders et al.l 1198911 . Variability 
in the accretion disk emission results in corresponding 
variations in the dust IR emission, but with a delay 
due to differing light travel times between the source, 
various points in the torus and the observer. The 
lags between the optical continuum light curve and 
the IR light curves can therefore be interpreted as 
measures of the distance from the source to the dust 
clouds that predominantly emit the 3.6 /im and 4.5 /rm 
radiation. Our results indicate these clouds are located 
at a distance w 1 light-month (« 0.03 pc) from the 
source of the ACN UV-optical continuum. However, 
the two Spitzer bands have significantly different lags, 
with the 4.5 ^m-optical lag being longer by 9.9 ± 3.9 
days. The lag between the 4.5 ^m and 3.6 /rm light 
curves is 13.9 ± 0.5 days and is consistent with this 


difference. This implies that the clouds producing the 
bulk of the 4.5 /rm emission are about 10 light-days 
( ~ 27%) further from the UV-optical continuum source. 

In most models, the innermost radius of the torus is 
taken to be the dust sublimation radius which, for a typ¬ 
ical IS M dust composit i on with silicate grains of average 
size, is (lBarvainislll987t iNenkova et al.ll2008bll 


where L^oi is the bolometric Inminosity of the ACN 
and Tsub is the dust sublimation temperature. 

However, many broad-line ACN exhibit a distinct 
near infrared ” bump”, peaking around 2 — 4^m, 
which has a black body t empe rat ure T > 1 000 K 
('e.g. lEdelson fc MalkanI 119861: iBarvainisI 119871: 
Rodriguez-Ardila fc Mazzalayl 12004 iRiffel et al l 


2009 


:igue 2 

mr 


This feature often dominates the NIR 
and it has been found that it cannot be reproduced by 
torus models alone in fits to the infrared spectral energy 
distribution (SED); instead, one must add a separate 
hot (T ~ 1400 K) black body component. The latter has 
been attri buted to hot pure graphite dust located withi n 
the torus (iMor et al.l 120091 : iMor fc Trakhtenbrod 1201 ill . 
and iMor fc Netzeir ( 2012l lliave modeled this component 
as dust embedded in the outermost BLR, between the 
sublimation radius for pure-graphite grains. 


Rd,c — 0.5 


Lboi \ 1800 a: V® 


lO"^® erg s 


-1 


\ Tsub J 


pc (3) 


and the torus inner radius as g iven by equatio n [2] The 
hot dust spectrum computed by IMor &: Netzed suggests 
that this hot graphite dust contributes significant Inmi¬ 
nosity at 3.6 /rm and 4.5 /xm. 

In order to estimate the sublimation radii given by 
equations [2] and in it is necessary to determine Lboi- How¬ 
ever, this is difficult to determine accurately for NCC 
6418, as the optical spectrum is dominated by the stel¬ 
lar contin uum and the ACN it self is evidently heavily 
reddened (|Remillard et al.iri993[l . 

The Sloan D igital Sky Survey (SDSS) optical spectrum 
of NCC 6418 (jAhn et al.ll2012h (Figure shows broad 
Ha and narrow lines of [OHIJA5007, Ha, [NH]A6548,6583 
and [SH]A6717,6731, but the continuum is dominated by 
an evolved stellar population. The fact that the broad 
H/3 line is not evident in the spectrum indicates a steep 
broad-line Balmer decrement a nd suggests classi f icatio n 
as a Seyfert Type 1.9 (Syl.9; lOsterbro"^ 119771 1198111 . 
However, NCC 6418 is unusual in that the narrow H/3 
emission is also very weak (in fact, this line appears in 
absorption) and the [OHIJA4959,5007 lines are much 
weaker relative to the stellar continuum than is typical 
in Seyferts, even Sy 1.9s. Interestingly, these lines are 
not obviousl y visible in the earlie r (1989) spectrum 
obtained bv iRemillard et id] ()1993[ 1. even though the 
broad Ha line is clearly much stronger relative to the 
narrow Ha and [Nil] lines than in the SDSS spectrum. 
Evidently, the strong stellar continuum, the foreground 
reddening and the variable broad emission lines make 
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Figure 3. Cross-correlation centroid distributions (CCCDs) for 3.6 /im versus 4.5 fim (top), 4.5 fim versus optical (middle), 3.6 fim versus 
optical (bottom). We have shifted the 3.6 fim versus 4.5 ^m CCCD by 37.2 days, approximately the time lag between of the 3.6 fim and 
optical light curves, since, in principle, we expect its peak to coincide with that of the 4.5 ^m versus optical CCCD. 
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Figure 4. SDSS DR9 optical spectrum of NGC 6418, 
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Figure 5. Fit to the Ha. broad emission line and blended nar¬ 
row lines in the SDSS DR9 optical spectrum of NGC 6418. The 
fitted gaussian profiles represent: broad Ha (red), [N IIJAA6548, 83 
(magenta), and [S IIJAA6717, 31 (cyan) and narrow Ha (green). 

the classification of this source somewhat ambiguous. 

To determine the bolometric luminosity of the AGN, 
we used the relationship established between the broad 
Ha luminosity and the bolometric AGN lumi- 

nosity {Lacn) a large sample of qu asars and Syl 
(|Richards et alJl^OOGt iStern fc Laoi]l201^ . 

Lioi = 130^1^ X (4) 

The flux in the broad Ha line was measured from the 
SDSS spectrum using gaussian profiles to fit and deblend 
the [N H], [S H] and Ha lines. In the fit, the wavelengths 
of the components representing the narrow lines were 
fixed at the values determined by the SDSS spectro scopic 
measurement pipeline (svecld: iBolton et al.fl201^ . The 
amplitudes and widths were free parameters, with the 
exception of [NH]A6548, which is constrained so as to 
preserve its fixed 1:3 intensity ratio with [NH]A6583. 
The variances provided by the SDSS spectroscopi c data 
reduction pipeline {spec2d\ iStoughton et al.l [200211 were 
used to assign weights to each data point; in addition, 
we ass igned a 10% syst ematic error to the derived 
fluxes (IBolton et al.l 120121) . The resulting fit is shown 
in Figure [S] and the parameters derived from the fit 
are summarized in Table SI The broad Ha component 
has a flux of (2563 ± 120) xlO”^"^ erg s“^cm“^. Using 
this Ha flux and assuming a distance of 122 Mpc 
(|Mould et al.ll2000l) . we calculate the observed Ha broad 
line luminosity to be L°^l, =(4.56 ± 0.85) x 10^° erg s“^. 

It is clear, however, that a large extinction correction 
needs to be applied in order to obtain the intrinsic Ha 
luminosity. From the SDSS spectrum we estimate a lower 
limit to the broad line Balmer de crement of Ha/H/3 > 6. 
We used the mean Ha/HP from iDong et al.l (1200511 and 
their expression to allow for reddening: 

\ogHSl = logL?,'^ + I.87(log(i7a/FI/3)-log(2.97)) (5) 

which yields a lower limit to the intrinsic broad Ha 
luminosity of Up/l > (1.70 ± 0.32) x lO^^ erg s'b 

With this lower limit, equation |4] yields a lower limit 
to the bolometric Inminosity of the AGN in NGC 6418 



V flux [mjy] 

Figure 6. Flux variation gradient diagram of NGC 6418 con¬ 
structed from observations made at the Fountainwood Observatory 
in Southwestern University. The data are represented by the black 
dots. The host contribution as indicated by the asterisk is 0.58 
mJy and 0.29 mJy for the V and B band, respectively. The dashed 
lines indicate the ra nge of host slopes determined in the optical by 
ISakata et al.l (120101) . The dot-dashed and solid lines indicate the 
least-square best fit to the range of the AGN slope. 

of Lbol > (2.2lti;29)xl0"^^ erg s“^. Using liTaspi et all 
(|2000ll ’s relation Lboi ^ 9ALa(5100) and assuming 
La(5500) ^ La(5100) we obtain a lower limit to the AGN 
V-band luminosity of Ly > (2.461^1)x 10^^ ergs“^. 
For comparison, we used the flux variation gradient 
(FVG) method d Choloniewskil 119 8^ ISakata et al.l [201(11 : 

iHaas et mi2011l : iPozo Nunez et al.ll2012L 1201411 to esti¬ 

mate the (constant) host galaxy contribution within our 
photometric aperture as illustrated in figured Using the 
B and V fluxes obtained from the SU observations (3.5” 
aperture), we And an AGN/Host ratio of 1.55, yielding 
an an estimate for the AGN contribution to the V-band 
luminosity of 1.54±0.53 x 10"^^ erg s“^ (as reddening cor¬ 
rections have not been applied to the B and V fluxes, this 
value should be regarded as a lower limit.) Thus, within 
the admittedly large uncertainties, the AGN V-band lu¬ 
minosity estimated from the FVG method is consistent 
with that determined from the Ha luminosity. 

Having determined the lower limit on the bolomet¬ 
ric luminosity of the AGN, we can determine the dust 
sublimation radii given by equations [5] and [31 For sil¬ 
icate dust with a sublimation temperature ~ 1500 K, 
we And Rd^Si > 60^2? x 10“^ pc ( 71 I 25 light days), 
whereas for pure graphite dust with sublimation temper¬ 
ature ~ 1800 K, we find Rd,c > 24^8^ x 10“^ pc (28tio 
light days). 

These sublimation radii bracket the radii derived from 
the lags at 3.6 ^m (Rr.s.e = 31.2l^ g x 10“^ pc) and 

4.5 pm (i?r,4.5 = 39 .5t2'6 ^ pc). As Rd^Si and 
Rd,c are lower limits, this suggests that the laulk of 
the 3.6 pm and 4.5 ^m emission comes from the region 
bounded by the graphite and silicate sublimation radii, 
respectively, and is conceivably emitted by the same 
graphite dust that is thought to be responsible for the 
NIR bump. As already noted, t he model graph i te dus t 
emission spectrum computed bv iMor &: NetzeJ (|2012f) . 
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V luminosity [erg/s] 


Figure 7. Reverbe r ation lag distanc e as a f unction of op t ical A GN luminosity. The data points are the if—band lag measurements of 
IKoshida et al.l Il2014lh FSuganuma et al] Il2006l h FClavel et al.l II1989I ) and the 3.6 /rm and 4.5 fim lag measurements of NGC 6418. The solid 
line represents the fit to the (t oc relationship as found by IlSuganuma et ani2006l 'l and defined as equation 3 in IIKishimoto et ahl 

1200711 . 


Table 4 

Emission line fit parameters 



line 


(A) 

fiux^ 

(ergs“^cm“^ x 10“^^) 

FWHM'^ 

(A) 

H 

Q broad 

6742±2 

2563±120 

156±4 

He 

narrow 

6753 

41±20 

5±1 


Nil 


6739 

90±25 

8±1 


Nil 


6775 

269±26 

8±1 


SII 


6912 

159±28 

10±1 


[5//] 


6926 

151±27 

9±1 


while peaking in the 2 — 3/im range, also emits strongly 
in the 3.6 - 4.5 fim range. Nevertheless, the longer lag 
exhibited by the 4.5 emission implies the presence 
of a temperature gradient in the emitting region. 


In K-band reverberation mapping studies of Seyfert 
1 galaxies it has been found that the reverberation ra¬ 
dius derived from the time lag is quite tightly correlated 
with where Lopt is the AGN optical luminosity 

l|Suganuma et al.ll2006t H^shida et ^1200911201411 . This 
is consistent with the R oc ® relat ion expected for 
dust i n radiative equilibrium. However. iKishimoto et ^ 
(|2007ll found that the AT—band reverberation radii are 
a factor 3 smaller than the sublimation radii as 
predicted by equation [2l One possible explanation is 
that the NIR dust emission is dominated by graphite 
grains; sublimation radii predicted by equation [3] are 
a factor ^ 3 smaller than the Silicate radii and thus 
much closer to the AT—band reverberation measurements 
(see Fig. ED. Several other explanations have been ad¬ 
vanced for the apparent discrepancy between the mea¬ 
sured dust radii and the sublimation radii predicted for 


the standard ISM dust composition. For example, the 
dust may include larger grains than the typical size 
(a PS 0.()5/xm) assumed in equa tion IIKishimoto et all 
1200711 . iKawaguchi fc Moril (1201011 investigated the effect 
of anisotropic illumination of the torus inner wall by the 
accretion disk, which permits a smaller torus inner radius 
close to the disk p l ane. A nother possibility, proposed by 
iPozo Nunez et all (|2014ll . is that the torus is very op¬ 
tically thick in the NIR so that only emission from the 
facing rim of the torus inner wall is seen, leading to a 
“foreshortened” lag. Modeling of the time-dependence 
of the optical -NIR spectral energy distribution (SED) of 
NGC 4151 bv iSchniille et ahl (j2(il3ll suggests that the in¬ 
nermost dust is well below the sublimation temperature. 
This implies that the dust is located beyond the subli¬ 
mation radius, suggesting anisotropic illumination or ge¬ 
ometr ical foreshortening, as envisaged iPozo Nunez et^ 

(IMl . 


In Figure El we plot reverberation radii versus V- 
band luminosity (ALa(V)) for both the 3.6 /rm and 
4.5 u m lags repo r ted h er e and K-band r e sults taken 
from iClavel et all (I1989D : iSuganuma et ahl (I2006D and 
IKoshida et ahl ( ^141) . For this purpose, we use the lower 
limit to the AGN V-band luminosity of NGC 6418 in¬ 
ferred from as described above. 

We also plot IKishimoto et al.l (|2007fl ’s fit to the K-band 
lag data points. 


Rt,K 


0.47 


6 A£y(V) 

1046 ergs“^ 



( 6 ) 


With the caveat that the NGC 6418 points represent 
lower limits in luminosity, it can be seen that the mid- 
IR reverberation radii are located above the trend de¬ 
fined by the AT-band lag times, as expected if the 3.6 
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fim and 4.5 emission is dominated by cooler dust lo¬ 
cated somewhat deeper in the torus. Equation |6] predicts 
Rt.k ^ 11-6 X 10“^ pc for NGC 6418, given our lower 
limit on the V luminosity, implying that i?r, 3.6 ^ 2.7i?r,K 
and i?r, 4.5 ^ 3.4i?.r,Kj respectively. 

For dust grains in radiative equilibrium, the radius at 
which grains have a temperature T is approximately. 


Rsub V Tsub J 

where Rsub is the sublimation radius and a ss 2 — 2.8 
depends on the dust composition. In combination with 
Wien’s Law, Equation [3 provides a rough estimate of the 
largest radius at which the dust contributes to the torus 
emission at a specific wavelength. For the typical ISM 
composition of Equation [3 (a = 2.6), we find Rs.q/Rk — 
3.6, Ra.b/Rk — 6.4 and R 4 . 5 /R 3.6 — 1.8. The values for 
Rs.g/Rk and Rss/Rk exceed the empirical upper limits 
determined from reverberation mapping, while the value 
of i? 4 . 5 /i? 3,6 agrees with the ratio of the reverberation 
lags (RrA-b/Rr, 3.6 = 1-3 ± 0.7). 

How ever, in clumpy torus models(e.g. iNenkova et al.l 
l2008allbl l. there is a wide range of dust temperature 
within a typical cloud, whic h therefore emits a broa d IR 
spectrum. In the models of INenkova et al.l (|2008bll . the 
bulk of the emission at A < 5/im emerges from clouds at 
no mor e than twice the inner radius (see INenkova et al.l 
l|2008bf) their Fig. 13). Thus, the relative sizes of the 
reverberation radii at 3.6 /rm, 4.5 fim and FT—band seem 
consistent with at least some clumpy tori models. 

It is well establi shed, mainly from Balmer line rever- 
beration mapping (|Bentz et al.ll2013l: IGreene et al.ll2010l 
and references therein) that the broad emission line re¬ 
gion follows a similar R ^ size-luminosity relation¬ 
ship. For a given AGN luminosity, the BLR reverbera¬ 
tion radius is a factor 4 — 5 small er than the K-band dust 
emission reverberati on radius (|Suganuma et al.l 120061 : 
iKoshida et al.l I2014I1 . as expected in the AGN uni- 
fication parad i gm. Interestingly, radii derived from 
IMor fc Netzed i|2012f) ’s SED fits suggest that the NIR 
emission component attributed to hot graphite dust 
clouds occupies a region intermedia te between the BLR 
and K-band reverberation radii fsee lKoshida et aT1l20l4 
Figure 13), consistent with the idea that this dust re¬ 
sides in the outer BLR clouds. In their analysis of mid- 
IR (12 ^m) interferometric observations, IBurtscher et al.l 
(1201311 find that although source sizes scale in a simi¬ 
lar way with Inminosity, there is a much larger scatter, 
with mid-IR source radii ranging from < 4 to 20 x i?r,K- 
A clearer picture of the structure of the AGN emission 
regions beyond the accretion disk is therefore beginning 
to emerge. Placing our results in this context, the re¬ 
verberation radii derived from the 3.6 fim and 4.5 fim 
light curves are consistent with the variable emission at 
these wavelengths arising in the inner clouds of the torus. 
However, we note as a caveat that NGC 6418 exhibits an 
atypical optical spectrum for a Seyfert 1, with a rela¬ 
tively strong, broad Ha line but with a steep Balmer 
decrement, relatively weak narrow lines (for instance, the 
equivalent width of [OHIJA5007 is only ^ 3A, that of 
narrow Ha ~ 0.5A) and with stellar emission dominat¬ 
ing the optical continuum. This indicates that the BLR 


and AGN UV-optical continuum are subject to heavy 
extinction along the line-of-sight, raising the possibility 
that the circum-nuclear dust distribution may be more 
quasi-spherical than toroidal in nature. 

5. SUMMARY 

We have presented initial results from the first year 
of a two-year campaign of IR (3.6 fim and 4.6 fim) 
and optical (B and V) monitoring of a sample of 12 
Seyfert 1 galaxies using the Spitzer Space Telescope 
supported by ground-based optical observations. In 
NGC 6418, we have found a lag between the mid-IR 
and optical light curves, with a time delay of 37 . 2^2 2 
days (31.2ti;9 x 10“^ pc) at 3.6 fim and days 

(39 .512 6 ^ 4.5 ^m, respectively. The 3.6 

fim emission leads the 4.5 fim emission by 13.9lgl days 
(11.7lgl ^ 4*^”^ P*^)- Tliese results indicate that the dust 
emitting the bulk of the 3.6 fim and 4.5 fim emission is 
located at a distance ss 1 light-month (r:^ 0.03 pc) from 
the source of the AGN UV-optical continuum. 

The nucleus of NGC6814 appears to be heavily red¬ 
dened, with a broad line Balmer decrement of Ha/H/3 
> 6 . For this reason, we can only determine a lower 
limit for the intrinsic luminosity of the AGN and hence 
lower limits on the dust sublimation radii. The reverber¬ 
ation radii are a factor ^ 2 smaller than the sublimation 
radius lower limit for silicate grains (sublimation tem¬ 
perature ~ 1500 K; Rd,Si > 6 OI 21 X 10“^ pc), but con¬ 
sistent with that for pure-graphite grains (sublimation 
temperature ~ 1800 K; i?d,c > 24jig^ x 10“^ pc). Rever¬ 
beration radii derived from K-band variability studies of 
other Seyferts are similarly a factor ~ 3 smaller than the 
silicate sublimation radius. It seems possible that some 
of the emission in the 3.6 - 4.5 fim range comes from hot 
graphite dust located within the region bounded by Rd,c 
and Rd,Si^ whose presence is suggested by SED model¬ 
fitting. 

The 3.6 and 4.5 fim reverberation radii fall above the 
extrapolated K-band size-luminosity relationship by fac¬ 
tors < 2.7 and < 3.4, respectively, while the 4.5 fim re¬ 
verberation radius is only 27% larger than the 3.6 fim ra¬ 
dius. This indicates a steeper temperature gradient than 
expected for optically thin dust in radiative equilibrium 
but is consistent with clumpy torus models, in which in¬ 
dividual optically thick clouds emit strongly over a broad 
wavelength range. 
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friend, colleague and mentor David Axon, who initiated 
this project and brought the collaboration together. We 
would like to thank Davide Lena for improvements to 
the manuscript and discussions on spectroscopic data re¬ 
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APPENDIX 

APPENDIX A - ENSEMBLE PHOTOMETRY 

We begin by defining a region within each image containing NGC 6418 and several nearby reference stars. Next, we 
determine a background value for this region by fitting a gaussian to the histogram of pixel values: the peak yields the 
background value and the width its uncertainty. Sources are detected usin g the STARS program of XVISTA, which 
employs an algorithm based on the FIND procedure within DAOPHOT (jStetsonI Il987f l. Candidate objects which 
survive cuts in several parameters such as full-width at half-maximum, sharpness and roundness, are selected for 
aperture photometry. We measure the brightness of each object using the PHOT program in XVISTA, which sums all 
counts within a circular aperture, including weighted contributions from pixels that lie partially inside the aperture. 
PHOT also measures the median pixel value within an annulus around each object to determine a local sky value and 
subtracts this from the object counts. Finally, the remaining object counts are converted to an instrumental magnitude. 

The se cond stage of the analysis subjects the measured instrumental magnitudes to inhomogeneous ensemble pho¬ 
tometry (|Honevcuttl 1199211 . Small differences in sky brightness, transparency, exposure time, and other factors can 
cause all objects in some particular exposure to appear slightly brighter or dimmer than average; ensemble photometry 
is designed to ident ify these systematic changes and remove their effects. 

iHonevciittl (jl992l l defines the equation of condition as 

m{e, s) = mO{s) + em{e), (Al) 

where m{e, s) is the instrumental magnitude of star s in exposure e and mO is the intrinsic instrumental magnitude 
of that star. The “exposure magnitude”, em, of an image accounts for variations in extinction, exposure time, 
background intensity and other effects that are common to all sources in an image. We note that even without the 
transparency issues that are typical of ground observations, the Spitzer IR data will have small variations due to 
changes in orientation and background illumination of the space telescope. The quantity that we want to minimize is 

ee ss 

P = ^ s) — mO(s) — em{e)]^w{e, s), (A2) 

e=l s=l 

where w{e,s) is the weight of each instrumental magnitude; we take its value to be CT(m(e, s))“^. This technique 
yields the best fit value of mO(s) for each source, assuming no intrinsic variability, and an empirical estimate of the 
uncertainty. In an ideal experiment, the uncertainty would be equal to that derived from the quadrature sum of the 
shot noise of the source, the sky noise and the detector read noise. This empirical estimate of the uncertainty is valid 
only for constant sources, such as the reference stars, but not for sources which vary intrinsically from one image to the 
next. The estimated uncertainties for the Spitzer Channel 1 data are shown in Figure [8] as a function of instrumental 
magnitude, with a quadratic fit to the reference stars in the field. 

APPENDIX B - CROSS-CORRELATION CODES RESULTS 

In an effort to give a comprehensive picture of the results obtained by different software packages used to determine 
the lag time between light curves we have included this appendix with table [5] of all results. The table contains the 
analyses of individual and combined optical datasets versus the infrared channels of the Spitzer Space Telescope. The 
first column indicates the Spitzer channel. The table has three sections, one for each of the software packages we used. 
Columns 2 through 6 are values obtained for our in-house cross-correlation package. Of those, columns 2-4 represent 
the difference from the median to the 25% value of the interquantile range (IQR), the median of the distribution 
and the difference from the median to the 75% value of the IQR, respectively. Columns 5 an d 6 are the mean an d 
the standard deviation. Column 7 is t he mean and the standard deviation for Peterson’s code (iPeterson et al.l[2003l . 
Columns 8-10 are Zu’s (iZu et al.ll20II[l corresponding to the low, mid and high values of the lag. 

APPENDIX C - CROSS-CORRELATION FUNCTION AND THE CROSS-CORRELATION CENTROID DISTRIBUTION 

Our simulations employ one thoudsand realizations of the light curves. Each synthetic light curve is generated 
by replacing each magnitude measurement with an artificial datum. This consists of the measured magnitude plus 
a random deviate drawn from a gaussian distribution with a mean of zero and standard deviation equal to the 
uncertainty in the measured value. We compute the CCFs and the corresponding weighted mean lags for each set of 
synthetic light curves to construct a distribution of the CCF centroids, the CCCD. 

The CCFs are often not symmetrical functions, and the skewness of these functions affects the calculation of their 
centroids. The question is - how to select the significant portion of each distribution, while discarding the uninteresting 
wings? Figure ini shows representative single realizations of the CCFs; it is obvious that the centroid of each CCF will 
depend on the range of data chosen for further calculation. In this work, we have adopted an algorithm that uses 
properties of each distribution itself to select the subset of measurements for the centroid calculation. First, we fit a 
cubic spline to the distribution in each realization, and compute the standard deviation, u, between the spline and 
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Ensemble Solution Uncertainty 



Figure 8. Spitzer chi (3.6 fim) uncertainties for the ensemble solution. The curve in red is the least-square fit to a quadratic function. 


Table 5 

Comparison of cross-correlation methods 


channel 

IQR 25% 

Vazquez 

median 

IQR 75% 

mean 

di std 

Peterson 
mean di std 

low 

Zu 

mid 

high 

dataset 

chi 

-5.1 

42.7 

5.3 

42.6 

± 

7.8 

42.3 ± 9.6 

-0.7 

40.9 

1.2 

FTN ISIS 

ch2 

-3.1 

50.4 

3.2 

50.5 

± 

7.1 

52.2 ± 13.5 

-6.5 

53.2 

5.2 

FTN ISIS 

chi 

-2.4 

35.0 

2.2 

35.1 

± 

3.7 

35.1 ± 4.2 

-0.8 

33.5 

6.9 

LT ISIS 

ch2 

-2.7 

47.3 

3.5 

47.8 

± 

4.0 

48.9 ± 4.4 

-1.2 

50.1 

0.6 

LT ISIS 

chi 

-5.5 

42.6 

5.4 

42.2 

± 

8.1 

47.5 ± 10.2 

-1.3 

39.6 

21.1 

FTN XVISTA 

ch2 

-3.6 

50.6 

3.0 

50.2 

± 

7.1 

60.2 ± 19.2 

-31.1 

69.0 

1.2 

FTN XVISTA 

chi 

-2.4 

27.1 

2.8 

27.5 

± 

4.1 

28.0 ± 6.0 

-0.7 

29.2 

0.8 

LT XVISTA 

ch2 

-2.6 

36.1 

3.3 

36.5 

± 

3.7 

35.2 ± 4.2 

-2.7 

33.2 

0.7 

LT XVISTA 

chi 

-2.4 

34.5 

2.1 

34.4 

± 

3.6 

34.5 ± 3.9 

-8.2 

35.6 

1.0 

LT -1- FTN XVISTA 

ch2 

-3.7 

44.6 

3.5 

44.5 

± 

4.2 

42.2 ± 4.6 

-2.7 

37.8 

6.8 

LT -1- FTN XVISTA 

chi 

-2.2 

37.2 

2.4 

37.2 

± 

3.3 

36.7 ± 3.4 

-6.5 

40.4 

0.7 

LT -1- FTN ISIS 

ch2 

-3.1 

47.1 

3.1 

47.3 

± 

4.6 

48.6 ± 3.7 

-4.7 

49.5 

1.2 

LT -1- FTN ISIS 

chl/ch2 

-0.1 

13.9 

0.5 

14.0 

± 

0.7 

14.6 ± 6.0 

-2.9 

13.2 

5.8 

CHI vs CH2 


the data. We adopt 2a as a measure of the dispersion within the CCF. We set a threshold in correlation which is the 
peak of the CCF minus this dispersion: K = peak — 2cr. All the CCF values greater than K are then used to calculate 
the centroid of that particular CCF. The fitted spline is shown together with the computed CCF(t) data points. We 
found that for optical vs 3.6 ^m, the top 24% of CCF data was used, for the optical vs. 4.5 ^m the top 23%, and for 
the 3.6 ^m vs 4.5 /im the top 6%. The threshold clearly is dependent on the noise characteristics of the underlying 
light curves which explains why the Spitzer light curves have a smaller data percentage used in the centroid calculation. 

After calculating the centroid of each realization of the CCF in this manner, we then combine all the centroids to 
create the cross-correlation centroid distribution (CCCD) for that pair of light curves. We choose the median value 
in the CCCD as the time lag between the two light curves. 
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CCF Spitzer CHI vs CH2 
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Figure 9. Sample realizations of the CCFs for each pair of light curves. The fitted line is a cubic spline. 





